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The radiative decays rj' (77) — > 77 are discussed. The modifications of the conventional PC AC formulae due to the 
gluoriic contribution to the flavour singlet axial anomaly are given. The decay constants satisfy a modified Dashen 
formula which generalises the Witten-Veneziano formula for the mass of the r\ . It is shown how the topological 
susceptibility in QCD with massive, dynamical quarks may be extracted from measurements of vj '(77) — > 77. 



1. The Result 

Radiative decays of the 77' are currently attract- 
ing renewed interest. (See j^j for references.) In 
this talk we summarise the special features that 
arise due to the gluonic contribution to the Ua (1) 
anomaly when PCAC methods (including chiral 
Lagrangians) are used in the flavour singlet chan- 
nel. In Q, we presented an analysis of 77' — > 77 
decay in the chiral limit of QCD, taking into 
account the gluonic anomaly and the associated 
anomalous scaling implied by the renormalisation 
group. Here, we summarise the results of a re- 
cent new analysis^) extending this to QCD with 
massive quarks, incorporating 77 — 77' mixing. In 
particular, we show how a combination of the ra- 
diative decay formula and a generalisation of the 
Witten-Veneziano mass formula for the 77' may 
be used to measure the gluon topological suscep- 
tibility x(0) in full QCD with massive quarks. 

Our main result is summarised in the formulae: 



f aa 5>j Q 77 + 2rifA g Gll 8 a0 = a", 



a 



(1) 



which describes the radiative decays, and 

f aa (m 2 Uf T ^ b = (2n f ) 2 A S a0 S bQ 
m u (uu) 
-2(/„„,tr T I m d (dd) 

m s (ss) 



(2) 
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which defines the decay constants appearing in 
(1) through a modification of Dashen's formula 
to include the gluon contribution to the J7a(1) 
anomaly. 

In these formulae, 77" denotes the neutral pseu- 
doscalars 7r°, 77, 77'. The (diagonal) mass matrix 
is (m 2 ) a i3 and 5^=77 is the appropriate coupling, 
defined as usual from the decay amplitude by 

(ll\v a ) = -igir-n e *pai3PiP2 eX (Pi) e (P?)- The 
constant cig m is the coefficient of the electromag- 
netic contribution to the axial current anomaly: 

a^Jfe = d Qc6 m c ^ + 277 / ,5 a og + < m ^^^F Alt/ (3) 



Here, J° 5 is the axial current, <j>% = qj 5 T a q is 
the quark pseudoscalar and Q = ^tvG^G^ is 
the gluon topological charge. 771° are the quark 
masses (see eq(ll)). a — 0,3,8 is the flavour 
index, T 3,8 are SU(3) generators and T° = 1. 
The d-symbols are defined by {T a ,T b } = d abc T c . 

The decay constants f aa in (1) are defined by 
the relation (2). In general they are not the cou- 
plings of the pseudoscalar mesons to the axial 
current |Q. In the flavour singlet sector, such a 
definition would give a RG non-invariant decay 
constant which would not coincide with the quan- 
tities arising in the correct decay formula (1). In 
contrast, all the quantities in the formulae (1),(2) 
are separately RG invariant 0,(3). The proof is 
not immediately obvious, and depends on the 
RGEs for the various Green functions and ver- 
tices defining the terms in (1),(2) being evaluated 
on-shell or at zero-momentum. 

In practice, since flavour SU(2) symmetry is al- 
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most exact, the relations for ir° decouple and are 
simply the standard ones with f 3v identified as 
/tt (see eqs(24),(27)). In the octet-singlet sector, 
however, there is mixing and the decay constants 
form a 2 x 2 matrix: 

aa ( f° V ' / 0T? "\ 

— ( J8))' J8v J (4) 

The four components are independent. In par- 
ticular, for broken 5t/(3), there is no reason to 
express f aa as a diagonal matrix times an orthog- 
onal 77 — 77' mixing matrix, which would give just 
three parameters. Several convenient parametri- 
sations may be made, e.g. involving two constants 
and two mixing angles, but this does not seem to 
reflect any special dynamics. 

The novelty of our results of course lies in the 
extra terms arising in (1) and (2) due to the glu- 
onic contribution to the !7a(1) anomaly. The co- 
efficient A is the non-perturbative number which 
specifies the topological susceptibility in full QCD 
with massive dynamical quarks. The topological 
susceptibility is defined as 

X(0) = J d 4 x i{0\T Q(x) Q(0)|0) (5) 

The anomalous chiral Ward identities determine 
its dependence on the quark masses and conden- 
sates up to an undetermined parameter, viz. 

X(0) = -A ( 1 - A Y —^7—) (6) 

Notice how this satisfies the well-known result 
that x(0) vanishes if any quark mass is set to 
zero. 

The modified flavour singlet Dashen formula is 
in fact a generalisation of the Witten-Veneziano 
mass formula for the 77'. Here, however, we do 
not impose the leading order in 1/N C approxi- 
mation that produces the Witten-Veneziano for- 
mula. Recall that this states 

6 

m v , + m v - 2m K = -^x(0)| YM (7) 

Jit 

To recover (7) from our result (see the first of 
eqs(9)) the condensate m s {ss) is replaced by the 
term proportional to f^m 2 K using a standard 



Dashen equation, and the singlet decay constants 
are set to y^2n~ff n . The identification of the large 
iV c limit of the coefficient A with the non-zero 
topological susceptibility of pure Yang- Mills the- 
ory follows from large N c counting rules and is 
explained in ref . 

The final element in (1) is the extra 'coupling' 
<7G 77 in the flavour singlet decay formula, which 
arises because even in the chiral limit the 77' is not 
a Goldstone boson because of the gluonic ?7a(1) 
anomaly. A priori, this is not a physical coupling, 
although (suitably normalised) it could be mod- 
elled as the coupling of the lightest predominantly 
glueball state mixing with 77'. However, this in- 
terpretation would probably stretch the basic dy- 
namical assumptions underlying (1) too far, and 
is not necessary either in deriving or interpret- 
ing the formula. In fact, the gcyy term arises 
simply because in addition to the electromagnetic 
anomaly the divergence of the axial current con- 
tains both the quark pseudoscalar 4% and the glu- 
onic anomaly Q. Diagonalising the propagator 
matrix for these operators isolates the 77 and 77' 
poles, whose couplings to 77 give the usual terms 
g vll and g^'^y However, the remaining operator 
(which we call G) also couples to 77 and therefore 
also contributes to the decay formula, whether or 
not we assume that its propagator is dominated 
by a ' glueball ' pole. 

Of course, the presence of the coupling <?g 77 in 
(1) appears to remove any predictivity from the 
if — » 77 decay formula. In a strict sense this is 
true, but we shall argue that it may nevertheless 
be a good dynamical approximation to assume 
(7G77 i s smau compared to <?,/ 77 . In this case, 
we can combine eqs (1) and (2) to give a mea- 
surement of the non-perturbative coefficient A in 
x(0)- To see this, assume m u = md — 0. The 
terms in (1), (2) involving the pion then decouple 
leaving the following five equations, in which we 
assume the physical quantities m v , m n ' , <?^ 77 and 
grfyy are all known and we neglect goy^- The 
decay equations are: 

/ I? S , r;7 7 + / 71 9r)'ll ~ "omj 

(8) 

where a° m = |iV c and a® m = jh^N c , and the 



3 



Dashen equations are: 



2. The Proof 



(/""') 2 m2, + (/°") 2 m2 = -4m s (ss) + 36 A 



(f 



8 n \ 



(/ 



m s (ss) 



(9) 



Clearly, the two purely octet formulae can be used 
to find / 17 and / 17 if both g m7 and g 7) ' 77 are 
known. The off-diagonal Dashen formula then 
expresses f 0r > in terms of /°'' . This leaves the 
two purely singlet formulae involving the still- 
undetermined decay constant / ^ , the topolog- 
ical susceptibility coefficient A, and the coupling 
gG-y-y- The advertised result follows immediately. 
If we neglect gG-y-y, we can find f 0r ' from the sin- 
glet decay formula and thus determine A from the 
remaining flavour singlet Dashen formula. This is 
the generalisation of the Witten-Veneziano for- 
mula. 

Without neglecting gc-yj, the five equations 
give a self-consistent description of the radiative 
decays, but are non-predictive. It is therefore im- 
portant to analyse more carefully whether it is 
really legitimate to neglect gG-y-y- The argument 
is based on the fact that gG-y-y is both OZI sup- 
pressed and renormalisation group (RG) invari- 
ant [0. Since violations of the OZI rule are as- 
sociated with the ?7a(1) anomaly, it is a plausi- 
ble conjecture that we can identify OZI-violating 
quantities by their dependence on the anomalous 
dimension associated with the non-trivial renor- 
malisation of J° 5 due to the anomaly. In this way, 
RG non-invariance can be used as a flag to indi- 
cate those quantities expected to show large OZI 
violations. If this conjecture is correct, then we 
would expect the OZI rule to be reasonably good 
for the RG invariant gG-y-y, which would therefore 
be suppressed relative to g^'-y-y- (An important 
exception is of course the rj mass itself, which 
although obviously RG invariant is not zero in 
the chiral limit as it would be in the OZI limit 
of QCD.) Notice that this conjecture has been 
applied already with some success to the 'proton 
spin' problem in polarised deep inelastic scatter- 
ing[Q. A related 1/N C discussion is given in [||. 



Consider first QCD by itself without the cou- 
pling to electromagnetism. The axial anomaly is 



d^J^ 5 = M ab r 5 +2n f Q6 a0 



(10) 



The notation is defined in rcf [Q] . The quark mass 
matrix is written as m a T a , so 







m u l + m 3 T 3 + m 8 T a (11 



The condensates are written as 
\ 

(dd) = -(</>°)l+2(0 3 )T 3 +2<0 8 )T 8 (12) 
(ss) J 6 

where (</> c ) is the VEV (qT c q). Then 




M, 



d a cb"m c 



dabM) (13) 

The anomalous chiral Ward identities, at zero 
momentum, for the two-point Green functions of 
these operators are|2j 

2n/(Q Q)S a0 + M ac {4>i Q) = 
2n f {Q 4>\)5 aQ + M ac {cf> c 5 4>\) + $ ah = 

(14) 

which imply 

M ac M bd {4>l 4) = -(M^) ab +(2n f ) 2 (Q Q)S a0 d b0 (U) 

We also need the result for the general form of 
the topological susceptibility (see eq(6)): 

X (0) = (Q Q}= 7 2 A A{Af ^-i ( 16 ) 

1 - (2n / )M(M$) 00 

Although the pseudoscalar operators 4>t an d Q 
indeed couple to the physical states r/ a = 7/ , 77, 7r°, 
it is more convenient to redefine linear combina- 
tions such that the resulting propagator matrix is 
diagonal and properly normalised. So we define 
operators r\ a and G such that 



(Q Q) {Q 4) 

($! Q) (05^) 

This is achieved by 

G = Q-{Q h))a b 4 
= Q + 2n / A$- b V| 



(G G) 






(r] a rf) 



(17) 



(18) 
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and 

V = J 



Taa^-l j.b 



b (19) 
With this choice, the (G G) propagator is 
(G G) = -A (20) 

and we impose the normalisation 
1 



(if rf) = 



k 2 



-6 a/3 



(21) 



This implies that the constants f aa in (19), which 
are simply the decay constants, must satisfy the 
(Dashen) identity 

f aa < P f Tpb =<M<&&»^*.a 

= -(M<S>) ab + (2n f ) 2 A6 a0 6 ba 

(22) 

The last line follows from the Ward identities (15) 
and (16) . In terms of these new operators, the 
anomaly equation (10) is: 



(23) 



3% a 5 = f aa m%if + 2n f GS a0 

Now recall how conventional PCAC is applied 
to the calculation of ir° — > 77. The pion decay 
constant is defined as the coupling of the pion to 
the axial current 

(0|J M 3 5 k)=^A (0\d^4 5 M = fnml (24) 

and satisfies the usual Dashen formula The next 
step is to define a 'phenomenological pion field' n 
by 

8»J% -» Uml-K (25) 

To include electromagnetism, the full anomaly 
equation is extended as in (3) to include the 
Fjiv contribution. Using (24) we have 

^<77l^ 5 |0) 
= Uml (77MO) + < m ^ (771^^10) 
= f n ml (tt vr)(77k) + a* m £( J7 \F^ F^\0) 

where (w w) is the pion propagator — l/(fc 2 — m 2 ). 
At zero momentum, the l.h.s. vanishes because of 
the explicit fc M factor and the absence of massless 
poles. We therefore find, 



7T 



(27) 



In the full theory including the flavour singlet 
sector and the gluonic anomaly, we find a similar 
result. The 'phenomenological fields' are defined 
by (23) where the decay constants satisfy the gen- 
eralised Dashen formula (22) . Notice, however, 
that they are not simply related to the couplings 
to the axial current as in (24) for the flavour non- 
singlet. We therefore find: 

l ^( 77 |J- 5 |0) = / aQ m2 /3 (7 7 |^|0) 

+2n / <77|G|0)<y a0 + < m ^( J7 \F^F^\0) 
= f aa <p(if ^(llW) 
+2n f (G G)( 7 7|G)<5 a0 + < m ^(77l^^|0) 

using the fact that the propagators are diagonal 
in the basis i] a , G. Using the explicit expressions 
(20), (21) for the Green functions, we find in this 
case: 



f aa 9v a n + 2n f A g Gll S a0 



a cm „ 
7T 



(29) 



where the extra coupling gc-y-y is defined through 
(28). This completes the derivation. It is evi- 
dently a straightforward generalisation of conven- 
tional PCAC with the necessary modification of 
the usual formulae to take account of the extra 
gluonic contribution to the axial anomaly in the 
flavour singlet channel, the key point being the 
identification of the operators i] a and G in (23). 

Finally, notice that these methods may equally 
be applied to other decays involving the rf such as 



'1 



U 7 where V is a light 1 meson, vf 



7T7T7, 



ip — ► 7/7, etc. 
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